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Abstract 

Resulting  from  a  Ph.D.  research  a  Vehicle  Simulation  Programme  (VSP)  is  proposed  and  continuously  developed.  It  allows  simulating 
the  behaviour  of  electric,  hybrid,  fuel  cell  and  internal  combustion  vehicles  while  driving  any  reference  cycle  [Simulation  software  for  com¬ 
parison  and  design  of  electric,  hybrid  electric  and  internal  combustion  vehicles  with  respect  to  energy,  emissions  and  performances,  Ph.D. 
Thesis,  Department  Electrical  Engineering,  Vrije  Universiteit  Brussel,  Belgium,  April  2000].  The  goal  of  the  simulation  programme  is  to 
study  power  flows  in  vehicle  drive  trains  and  the  corresponding  component  losses,  as  well  as  to  compare  different  drive  train  topologies.  This 
comparison  can  be  realised  for  energy  consumption  and  emissions  as  well  as  for  performances  (acceleration,  range,  maximum  slope,  etc.). 

The  software  package  and  its  validation  are  described  in  [J.  Automot.  Eng.,  SAE IEE  215  (9)  (2001)  1043L].  Different  hybrid  and  electric 
drive  trains  are  implemented  in  the  software  [Views  on  hybrid  drive  train  power  management  strategies,  in:  Proceedings  of  the  EVS-17, 
Montreal,  Canada,  October  2000].  The  models  used  for  the  energy  sources  like  fuel  cells,  batteries,  ultracapacitors,  flywheels  and  engine- 
generator  units  will  be  discussed  in  this  paper  in  three  stages:  first  their  functionality  and  characteristics  are  described,  next  the  way  these 
characteristics  can  be  implemented  in  a  simulation  model  will  be  explained  and  finally  some  calculation  results  will  illustrate  the  approach. 

This  paper  is  aimed  to  give  an  overview  of  simulation  models  of  energy  sources  for  battery,  hybrid  and  fuel  cell  electric  vehicles. 
Innovative  is  the  extreme  modularity  and  exchangeability  of  different  components  functioning  as  energy  sources.  The  unique  iteration 
algorithm  of  the  simulation  programme  allows  to  accurately  simulate  drive  train  maximum  performances  as  well  as  all  kind  of  power 
management  strategies  in  different  types  of  hybrid  drive  trains  [IEEE  Trans.  Veh.  Technol.,  submitted  for  publication]. 

©  2003  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  pollution  caused  by  transport  is  a  heavy  burden,  es¬ 
pecially  in  urban  areas  [5].  The  introduction  of  clean  vehi¬ 
cles,  like  electrically  driven  vehicles,  would  be  an  interesting 
move  in  the  direction  of  a  significant  reduction  in  harmful 
exhaust  gases,  with  a  view  to  a  sustainable  transport  policy 
[6] .  Since  several  years  the  family  of  the  battery  electric  ve¬ 
hicles  is  extended  with  the  hybrid  electric  and  recently  also 
with  the  fuel  cell  electric  vehicles.  These  vehicle  technolo¬ 
gies  use  different  energy  sources,  which  will  be  discussed 
in  this  paper. 
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The  energy  sources  that  will  be  considered  in  this 
manuscript  are  mainly  used  in  “series”  hybrid  vehicles.  In 
series  hybrid  drives  there  are  no  mechanical  connections 
between  the  internal  combustion  engine  and  the  wheels  [7]. 
All  thermal  energy  is  converted  first  into  mechanical  energy 
in  a  thermal  engine  and  subsequently  into  electrical  energy 
by  a  generator  driven  by  the  thermal  engine.  Additionally, 
there  is  an  electric  traction  motor  to  drive  the  wheels.  Hence 
a  decoupling  of  energy  source  operation  from  the  required 
traction  power  is  possible.  In  most  cases  the  energy  source, 
also  called  auxiliary  power  unit  (APU),  will  act  as  base 
power  unit  delivering  power  to  the  battery  or  directly  to 
the  electric  traction  motor.  While  driving  the  battery  acts 
as  peak  power  unit  or  energy  buffer.  The  series  hybrid  has 
the  advantage  of  operating  a  thermal  engine  in  a  selected 
optimal  operating  field,  for  instance,  with  low  specific  fuel 
consumption  in  the  torque-speed  operating  area  [8]. 

The  fuel  cell  hybrid  structure  is  a  series  structure  in  which 
the  engine  and  the  generator  are  replaced  by  a  fuel  cell 
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Fig.  2.  Vehicle  simulation  programme. 
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system  producing  electric  energy  starting  from  stored  hy¬ 
drogen  or  from  a  fuel  tank  feeding  a  reformer  to  produce 
hydrogen  [9].  The  excess  of  electricity  produced  by  the  fuel 
cell  can  be  stored  in  a  buffer  battery.  When  the  battery  is 
left  out  one  has  no  longer  a  hybrid  vehicle  but  a  fuel  cell 
electric  vehicle.  In  this  case,  the  fuel  cell  has  to  demonstrate 
enough  dynamics  to  meet  variable  power  demand. 

Fig.  1  illustrates  a  simulation  model  of  a  series  hybrid  with 
peak  power  unit,  where  the  generator  group  provides  a  con¬ 
stant  average  power,  the  battery  works  as  an  energy  buffer 
and  the  flywheel  (or  ultracapacitor)  caters  for  the  brake  and 
acceleration  power  peaks.  In  this  configuration  the  battery 
can  be  smaller  or  it  can  possibly  be  omitted.  The  efficiency 
and  reliability  of  the  flywheel  or  ultracapacitor  are  an  im¬ 
portant  factor  for  the  development  of  these  drive  trains. 

This  model  is  part  of  a  vehicle  simulation  programme 
(VSP)  [2-4]  whose  front  panel  is  illustrated  by  Fig.  2  and  for 
which  it  is  necessary  to  introduce  models  for  all  components 
or  subsystems  of  each  simulated  vehicle  and  particularly  of 
the  energy  and  power  sources.  The  basic  modelling  strategy 
used  in  VSP  is  the  well-tried  and  trusted  method  of  dividing 
the  drive  cycle  into  a  number  of  time  steps  and  calculating 
the  characteristics  of  the  vehicle  at  the  end  of  each  time  in¬ 
terval,  which  is  called  the  longitudinal  dynamics  simulation 
[10,11]. 

This  paper  will  focus  on  the  description  of  different  energy 
or  power  sources  as  well  as  on  their  simulation  models  and 
calculation  results. 

2.  Fuel  cells 

2.7.  Description  of  the  fuel  cell 

The  fuel  cell  electrochemically  combines  the  oxygen  of 
air  with  a  hydrogenated  fuel  and  converts  them  into  water 
and  other  elements  with  a  production  of  electricity.  The  con¬ 
version  is  similar  to  that  of  a  conventional  battery,  except 
that  the  reductant  and  oxidant  are  continuously  supplied  to 
the  cell  instead  of  being  contained  in  the  cell.  In  addition, 
fuel  cells  are  ‘recharged’  by  filling  up  the  fuel  supply. 

When  the  fuel  cell  uses  hydrogen  as  a  fuel,  this  hydrogen 
can  either  be  stored  as  such  in  the  vehicle,  or  be  generated  in 
the  vehicle  by  means  of  a  reformer.  Fuels  such  as  methanol, 
natural  gas,  synthesis  gas  or  even  conventional  hydrocarbon 
fuels  can  be  used  by  employing  a  chemical  reactor  prior  to 
entering  the  fuel  cell.  The  fuel  is  converted  into  a  hydro¬ 
gen  rich  synthesis  gas  via  a  steam  reformer,  auto-thermal 
reformer  or  a  partial  oxidising  burner.  This  gas  is  then  puri¬ 
fied  in  a  gas  purification  stage  to  clean  the  hydrogen  and  to 
scrub  CO  and  CO2,  which  are  emitted  from  the  process  [12]. 

Hydrogen  can  be  stored  as  a  compressed  gas  in  vessels, 
as  a  cryogenic  liquid  at  a  temperature  of  —253  °C  or  as  hy¬ 
drogen  atoms  reversibly  adsorbed  in  metal  hydrides,  which 
are  obtained  by  applying  hydrogen  gas  under  pressure  on 
suitable  metal  alloys.  The  low  temperature  and  pressure,  at 


which  the  reaction  occurs  and  the  stability  of  the  hydrides, 
preventing  leakage  flows  in  case  of  damage,  make  this  form 
of  storage  particularly  attractive  for  road  transportation  ap¬ 
plications.  Compared  to  lead  acid  batteries  this  storage  con¬ 
tains  up  to  10  times  more  energy  per  unit  of  weight  [13]. 

Fuel  cells  have  the  potential  to  be  more  efficient  than  a 
gasoline  or  diesel  engine.  Fuel  cells  are  at  the  opposite  of 
most  electrical  devices  since  peak  efficiency  occurs  at  mini¬ 
mum  load:  efficiency  is  60%  at  no  load  and  40%  at  full  load 
at  the  current  state  of  development,  the  theoretical  maxi¬ 
mum  being  83%  [14].  The  fuel  cell  system  has  relatively  few 
moving  parts  (reactant  pumps  and  valves),  which  produce 
very  little  friction.  The  fuel  is  oxidised  at  low  temperature 
instead  of  being  combusted  at  very  high  temperatures  as  in 
an  internal  combustion  engine.  Low  emissions  are  a  result 
of  low-temperature  oxidation  and  appropriate  fuel  selection 
when  a  reformer  is  used  with  hydrocarbon  fuels.  If  methanol 
is  used  to  generate  hydrogen  with  a  reformer,  25  vol.%  of 
the  gas  output  is  CO2  [15].  In  the  case  of  hydrogen  (the 
ideal  fuel),  the  only  by-product  is  water. 

Fuel  cell  systems  (e.g.  PEMFC)  are  under  development 
and  still  relatively  large  and  heavy  (also  due  to  the  large 
volume  of  ancillary  equipment)  compared  to  conventional 
engines.  Cost  is  still  very  high,  mainly  because  of  the  use  of 
expensive  materials  and  the  low  volume  production.  As  man¬ 
ufacturing  issues  are  tackled,  lower  cost  materials  should 
considerably  bring  the  cost  down. 

The  fuel  distribution  infrastructure  for  hydrogen  is  very 
limited.  A  reformer  based  vehicle  system  has  the  advantage 
of  using  existing  hydrocarbon  fuels  that  are  more  readily 
available  than  hydrogen: 

•  Methanol  is  a  simple  molecule,  making  it  relatively  easy 
to  reform  and  it  is  liquid,  so  it  fits  more  easily  into  the 
existing  motor  fuel  distribution  infrastructure.  Concerns 
are  raised,  however,  on  its  toxicity;  the  fact  that  it  burns 
without  readily  visible  flames  may  also  create  a  potential 
hazard. 

•  When  using  gasoline  reformers  the  existing  fuel  infras¬ 
tructure  can  be  used,  but  their  additional  emissions  are 
still  being  examined. 

•  Hydrogen  can  also  be  produced  out  of  natural  gas  at  re¬ 
fuelling  stations  [16]  or  by  electrolysis  of  water. 

The  fuel  cell’s  working  principle  is  very  simple,  how¬ 
ever,  the  electrochemical  reactions  do  not  readily  take  place, 
and  unless  special  materials  are  used  for  the  electrodes  and 
the  electrolyte,  the  current  produced  per  square  centimetre 
is  small  and  the  electrical  power  losses  in  the  electrolyte 
are  large.  To  overcome  these  problems,  different  types  of 
fuel  cells  have  been  developed.  The  different  varieties  are 
distinguished  by  the  electrolyte  used  and  have  each  a  par¬ 
ticular  construction  of  the  electrodes.  However,  in  all  types 
there  are  separate  reactions  at  the  anode  and  the  cathode, 
and  charged  ions  move  through  the  electrolyte,  while  elec¬ 
trons  move  through  the  external  load  circuit.  Another  com¬ 
mon  feature  is  that  the  electrodes  should  be  porous,  because 
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Table  1 


Summary  of  fuel  cell  types  [18] 


FC  type 

Electrolyte 

Operating 
temperature  (°C) 

Current 

density 

Need  for  fuel 

processor 

Stage  of  development 

Efficiency 

prospects 

Cost 

prospects 

PEMFC 

PEM 

70-90 

High 

Yes 

Prototype  demonstrated  on 
board  of  vehicle  fleet 

Good 

Good 

AFC 

Alkaline 

70-100 

High 

Yes 

Space  application 

Good  (70%) 

Good 

PAFC 

Phosphoric  acid 

200-220 

Moderate 

Yes 

Early  commercial  (stationary) 
applications 

Good  (40%) 

Fair 

MCFC 

Molten  carbonate 

600-650 

High 

Yes 

Field  demonstration,  high 
power  applications 

Good  (60%) 

Fair-good 

SOFC 

Solid  oxide 

800-1000 

Moderate 

No 

Prototype 

Good  (60%) 

Fair-good 

DMFC 

PEM 

70-90 

Moderate 

No 

Research 

Poor  (40%) 

Poor-fair 

the  gases  have  to  be  in  contact  with  the  electrode  and  the 
electrolyte  at  the  same  time  [17].  Table  1  summarises  the 
features  of  the  different  fuel  cell  types. 

2.2.  Fuel  cell  model 

In  the  previous  section  the  basic  principles  of  fuel  cell 
are  described.  However,  how  can  they  be  expressed  in  a 
simulation  model?  This  will  be  explained  hereunder. 

The  fuel  cell  model  considered  is  a  typical  black  box  ap¬ 
proach.  This  makes  it  possible  to  develop  in  the  first  in¬ 
stance  one  simplified  model  containing  the  fuel  cell  stack, 
the  DC-DC  convertor,  the  reformer  and  the  auxiliary  system 
[19]. 

Up-to-now  temperature  and  pressure  effects  are  not  in¬ 
cluded  in  the  model,  due  to  lack  of  sufficient  data  avail¬ 
able.  A  higher  stack  temperature  or  lower  pressure  will 
lead  to  a  lower  voltage  level.  The  start-up  process  of  some 
fuel  cells  can  be  very  long  and  can  create  considerable  ad¬ 
ditional  energy  consumption  and  emissions.  The  model  is 
based  on  experimental  data  of  an  alkaline  fuel  cell  (AFC), 
as  well  as  of  a  proton  exchange  membrane  (PEM)  fuel  cell 

[15]. 

The  DC-DC  convertor  adapts  the  fuel  cell  voltage  level 
to  the  DC-bus  voltage  or  battery  voltage.  At  the  moment  it  is 
modelled  by  a  simple  efficiency  value  (see  (2.1)).  The  power 
required  for  the  fuel  cell  auxiliaries,  Paux  (pump,  compres¬ 
sor,  controller)  is  calculated  based  on  experimental  data  in 
function  of  the  required  fuel  cell  power,  Pfc  (see  Eq.  (2.2)) 


[15]: 

Pdc 

Pic  = 

Oconv 

(2.1) 

Aux  —  O  +  bP  fc  +  cP 

(2.2) 

Aef  =  Yref  Pfc 

(2.3) 

P stack  =  Ac  T  Aux  T  Aef 

(2.4) 

where  Pdc  is  the  DC/DC  convertor  power  to  DC-bus  (W); 
Tjconv  the  DC/DC  convertor  efficiency  (%);  Pfc  the  fuel  cell 
power  (W);  a ,  b ,  c  the  constant  empirical  values;  Paux  the 


auxiliary  power  (W);  vref  the  reformer  coefficient;  Aef  the 
reformer  power  (W)  and  Atack  the  fuel  cell  stack  power  (W). 

The  power  required  for  the  reformer,  Aef>  is  also 
taken  into  account  as  proportional  to  the  delivered  power 
(Eq.  (2.3)).  The  sum  of  these  power  components  results  in 
the  total  required  stack  power,  Atack  (Eq.  (2.4)).  If  this  stack 
power  exceeds  the  maximum  allowed  stack  power,  Anax, 
the  iteration  process  [20]  of  the  simulation  programme  will 
decrease  the  fuel  cell  power,  Pfc. 

The  fuel  cell  stack  consists  of  different  modules  connected 
in  series.  Each  module  is  a  combination  of  a  number  of 
individual  cells  in  series  and  different  rows  of  series-cells  are 
connected  in  parallel.  The  required  power  for  each  cell,  Aelk 
can  be  calculated  by  dividing  the  stack  power  by  the  total 
number  of  cells.  This  is  an  approximation,  however,  due  to 
the  possible  unbalance  between  cells.  Based  on  experimental 
data  the  cell  current,  7ceii,  can  be  computed  as  function  of 
this  cell  power  [15]. 

In  the  considered  electrical  circuit,  the  fuel  cell  is  mainly 
characterised  by  its  internal  resistance.  The  dynamic  be¬ 
haviour  of  the  fuel  cell  has  also  to  be  taken  into  account: 
the  maximum  rate  of  change  of  the  current  is  dependent  on 
the  cell  surface.  The  simulation  model  shall  check  the  max¬ 
imum  current  variation  and  if  necessary  it  will  reduce  the 
required  fuel  cell  power. 

The  fuel  consumption  per  cell  (fee)  is  calculated  as  func¬ 
tion  of  the  current  law  of  Faraday,  which  gives  the  relation 
between  the  reactant  consumption  and  the  produced  current: 
1  g  equivalent  weight  of  matter  is  electrochemically  altered 
at  each  electrode  for  96,489  C  of  electricity  passed  through 
the  electrolyte.  Two  gram  of  hydrogen  (H2)  produces 
192,978 Ampere-seconds  (As)  or  53.61  Ampere-hours 
(Ah).  The  elimination  (purge)  of  impurities  is  taken  into 
account  with  a  correction  factor  (e.g.  1.012  in  Eq.  (2.5)) 
[15].  The  fuel  consumption  corresponds  with  a  certain  fuel 
power  content  as  function  of  the  density,  D  (kg/m3),  and 
the  specific  energy,  £h?  (Eq.  (2.7)): 


f  /cell  22.4  (l/mol) 

fee  =  1.012 - 

2.96489  (As/mol) 


(2.5) 


Cons  =  Nrceiis  fee 


(2.6) 
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Pfuei  =  Cons  Eu2  D  X  If)6  (2.7) 

where  /ceii  is  the  cell  current  (A);  fee  the  fuel  consumption 
per  cell  (1/s);  Cons  the  fuel  consumption  of  the  whole  stack 
(1/s);  Nrcens  the  number  of  cells;  D  the  fuel  density  (kg/m3); 
the  specific  energy  (MJ/kg)  and  Rfuei  the  fuel  power 

(W). 

The  emissions  are  calculated  proportional  to  the  fuel  con¬ 
sumption  [21]  of  the  related  reformer.  This  is  a  very  sim¬ 
plified  approach  that  only  gives  an  indication  of  the  order 
of  magnitude  of  the  emissions.  Additional  data  has  to  be 
available  for  the  extension  of  the  model.  This  model  is  an 
unpretentious  model,  but  due  to  the  black  box  approach  it 
can  easily  be  completed  as  soon  as  more  experimental  data 
are  available,  for  example,  out  of  the  ELEDRIVE  thematic 
network  [22]. 

2.3.  Calculation  results  of  the  fuel  cell  model 

Fig.  3  illustrates  the  simulation  result  of  a  fuel  cell  used 
in  a  2.7 1  hybrid  electric  van  driving  an  ECE-15  speed  cycle. 
The  input  parameters  are  the  requested  fuel  cell  power  as 
well  as  the  composition  of  the  fuel  cell  stack  (number  of 
individual  cells).  The  graphs  in  Fig.  3  illustrate  the  voltage, 
current  and  power  delivered  by  one  cell  of  the  fuel  cell 
stack. 

In  this  example  the  simulated  power  management  strategy 
of  the  hybrid  drive  train  consists  in  continuously  charging 
the  battery  at  5  kW,  while  the  fuel  cell  delivers  this  charging 
power  as  well  as  the  power  requested  by  the  electric  traction 
motor  for  driving.  This  somewhat  exotic  power  management 
was  selected  to  better  illustrate  the  dynamics  in  the  power 
delivered  by  the  fuel  cell.  Another  strategy  could  consist  in 
delivering  constant  power  by  the  fuel  cell,  the  battery  being 
in  charge  of  all  power  dynamics. 


3.  Battery 

3.1.  Battery  description 

Batteries  are  characterised  by  their  life  cycle,  energy  and 
power  density  and  energy  efficiency.  The  life  cycle  repre¬ 
sents  the  number  of  charging  and  discharging  cycles  possible 
before  it  loses  its  ability  to  hold  a  useful  charge  (typically 
when  the  available  capacity  drops  under  80%  of  the  initial 
capacity).  Life  cycle  typically  depends  on  the  depth  of  dis¬ 
charge.  When  charging  and  discharging  a  battery  not  all  en¬ 
ergy,  delivered  to  the  battery,  will  be  available  due  to  battery 
losses,  which  are  characterised  by  the  battery  efficiency.  The 
specific  energy  and  power  describe,  respectively,  the  energy 
content  (determining  the  vehicle  range)  and  the  maximum 
power  (determining  the  vehicle  acceleration  performance)  in 
function  of  the  weight  of  the  battery.  A  battery  can  be  opti¬ 
mised  to  have  a  high  energy  content  or  it  can  be  optimised  to 
have  a  high  power  capability.  The  first  optimisation  is  impor¬ 
tant  for  battery  electric  vehicles,  while  the  second  is  mainly 
required  for  hybrid  drive  trains.  Hybrid  drive  trains  lead  also 
to  the  consideration  of  a  mixed  energy-power  capability. 

A  specific  energy  of  at  least  50Wh/kg  and  a  specific 
power  ranging  from  lOOW/kg,  continuous,  up  to  200  or 
300W/kg  intermittent  pulse  power  is  a  must  for  an  good 
electric  vehicle  design.  Ideally,  150  W  h/kg  and/or  500  W/kg 
up  to  1000  W/kg  (in  the  case  of  super-capacitor)  would  be 
better,  especially  for  hybrid  vehicles  [9].  The  US  Advanced 
Battery  Consortium  (US ABC)  midterm  goals  are  150  and 
135  W  h/kg  and  long-term  goals  400  and  300  W  h/kg  for, 
respectively,  battery  power  and  energy  density. 

As  is  usual  in  electric  circuit  theory  the  maximum  power 
is  defined  for  the  system  delivered  rate  current  at  a  terminal 
voltage  that  equals  50%  the  initial  open  circuit  voltage  [23]. 


Fig.  3.  Simulation  result  of  a  fuel  cell  model. 
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In  the  IEC  standard  [24]  the  maximum  deliverable  battery 
power  is  defined  as  the  power  corresponding  to  rated  current 
and  the  battery  terminal  voltage  would  fall  at  2/3  of  the  open 
circuit  voltage. 

A  battery  for  an  electric  vehicle  consists  of  a  combina¬ 
tion  in  series,  and  possibly  parallel,  of  individual  battery 
cells.  Minor  differences  in  cell  characteristics  can  be  present. 
Consequently,  a  battery  management  system  controlling  the 
state  of  each  cell  is  a  key  point  for  battery  reliability  and 
life  duration. 


3.2.  Battery  model 


The  model  described  is  the  one  of  a  lead  acid  battery.  Sim¬ 
ilar  models  for  other  batteries  are  possible.  The  model  tracks 
the  state  of  charge  (SoC)  according  to  the  current  drawn 
from  the  battery.  The  SoC  is  the  relative  remaining  capac¬ 
ity,  C/,  of  the  battery,  expressed  in  %.  Discharging  during  a 
period  Ts  with  a  current  4at  leads  to  a  SoC  equalling: 


SoC  =  1 


A>  at  Ts 

C/3600 


(3.1) 


The  C5 -capacity  is  the  number  of  Ampere-hours  (Ah), 
which  can  be  drawn  out  of  the  battery  during  5h  with  a 
current  equal  to  C5  divided  by  5  (i.e.  one-fifth  of  the  rated 
Ampere-hours  in  5  h  discharge). 

The  effect  of  ambient  temperature,  t®  (in  °C)  on  the  ca¬ 
pacity  is  taken  into  account  with  the  help  of  an  empirical  for¬ 
mula.  An  example  for  a  lead  acid  battery  is  give  by  Eq.  (3.2) 
[25].  The  influence  of  age,  number  of  discharge/charge  cy¬ 
cles  are  yet  not  taken  into  account.  The  C5 -capacity  has  to 
be  multiplied  by  the  number  of  series-cells  in  parallel,  Nrp, 
to  get  the  total  battery  capacity,  Cto t' 


NrpC5 
1  +  0.006(30 


(3.2) 


The  available  capacity  of  the  battery  changes  as  function 
of  the  discharge  (or  charge)  current.  The  implemented  model 
tracks  the  SoC  according  to  the  current  drawn  from  the  bat¬ 
tery.  On  the  bases  of  the  actual  battery  current  /bat  and  the 
time  increment  Ts ,  the  variation  of  SoC,  A  SoC  is  calculated. 
It  is  primarily  based  on  the  empirical  law  of  Peukert  [26,27]. 

Peukert’ s  equation  is  a  formula  that  shows  how  the 
available  capacity  of  a  battery  changes  according  to  the 
rate  of  discharge.  The  capacity  of  a  battery  is  expressed  in 
Ampere-hours  (Ah),  but  it  turns  out  that  the  simple  for¬ 
mula  of  current  times  hours  does  not  accurately  represent 
the  situation.  Peukert  found  that  Eq.  (3.5)  fits  better  the 
observed  behaviour  of  batteries: 

C  =  Ik  Time  (3.5) 


The  equation  captures  the  fact  that  at  higher  currents,  there 
is  less  available  energy  in  the  battery.  C  is  the  theoretical 
capacity  of  the  battery  and  /  the  current,  Time  the  total  charge 
time,  and  k  the  Peukert  number,  a  constant  for  a  given  battery. 
A  value  close  to  1  indicates  that  the  battery  performs  well. 
The  higher  k ,  the  more  capacity  is  lost  when  the  battery  is 
discharged  at  high  currents.  The  Peukert  number  of  a  battery 
is  determined  empirically  and  lies  normally  between  1 .2  and 
1 .4  for  lead  acid  batteries.  If  k  would  be  selected  equal  to  1 ,  a 
linear  relation  between  SoC  and  capacity  would  be  assumed, 
as  this  is  practically  the  case  for  NaNiCh  batteries.  Fig.  4 
represents  this  equation  for  different  values  of  k.  It  is  quite 
remarkable  that  discharging  at  0.1  times  the  nominal  current 
increases  the  capacity  with  a  factor  2  when  the  battery  has 
a  Peukert  constant  k  equals  1.3. 

When  the  5-h  rate  capacity,  C5,  is  given  the  Peukert  ex¬ 
pression  can  be  used  to  calculate  the  capacity  for  any  given 
current: 


C/  =  C5 


k- 1 


(3.6) 


Within  this  model  a  subprogramme  characterising  the 
cell-module  is  used  to  calculate  the  open  circuit  voltage 
of  the  battery  pack,  I/o,  and  the  internal  resistance, 
both  dependent  on  the  current  SoC,  the  temperature  and  on 
discharge  or  charge  condition  (Eqs.  (3.3)  and  (3.4)).  The  de¬ 
pendence  of  the  internal  resistance  of  a  lead  acid  battery  on 
the  SoC,  for  example,  is  due  to  the  variation  of  electrolyte 
conductivity  and  of  the  bad  conducting  lead  sulphate  which 
deposed  on  the  electrodes  during  discharge.  The  dependency 
on  the  internal  resistance  on  the  SoC  is  bound  to  the  battery 
type.  These  figures  can  be  measured  for  charging  as  well 
as  for  discharging.  Notice  that  the  model  does  not  take  into 
account  the  possible  unbalance  between  cells  or  monoblocs. 

u0  =  Nrs/(SoC,  t° )  (3.3) 


Nrs 

Nr? 


/(SoC,  f°) 


(3.4) 


where  Nrs  is  the  number  of  cells  in  series;  Nrp  the  number 
of  series-cells  in  parallel;  t°  the  cell  temperature  (°C);  Uo 
the  open-loop  voltage  (V)  and  R[  the  internal  resistance  (£2). 


The  substitution  of  (3.6)  in  (3.1)  gives: 


SoC  =  1 


/bat  Ts  /  /bat  \ 

C~5  3600  [iTJ 


(3.7) 


This  law  is  transformed  into  (3.8),  which  expresses  the  SoC 
increment,  A  SoC.  This  equation  assumes  that  the  current 
/bat  remains  constant  during  one  step  Ts  of  the  simulated 
speed  cycle.  In  this  equation,  Ctot  (Eq.  (3.2)),  taking  into 
account  the  temperature  effect,  replaces  C5 : 


ASoC  = 


T,  14a/ 
3600  (Cm/S^-'Cm 


(3.8) 


To  be  able  to  implement  this  formula  the  knowledge  of 
the  battery  current  is  required.  This  value  can  be  calculated 
based  on  the  fundamental  electric  equations  (3.9)  and  (3.10), 
which  describe  the  simple  battery  model  of  Fig.  5.  Trans¬ 
forming  these  formulas  gives  (3.11),  which  defines  the  bat¬ 
tery  current  /bat  as  function  of  the  required  battery  power 
Ebao  the  internal  resistance  R[  and  the  open-loop  voltage  Uq: 


U bat  —  Uq  R\  /bat 


(3.9) 
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Fig.  4.  The  battery  capacity  vs.  discharge  current. 


^bat  —  ^bat-^bat 


^bat  — 


Uo  -  yjul  -  4/W/?i 


2R{ 


(3.10) 


(3.11) 


Dynamic  modelling  of  a  battery  is  extremely  difficult.  Un¬ 
derstanding  what  is  going  on  in  batteries  requires  several 
years  of  thorough  research  effort.  Finding  the  right  equation 
to  describe  these  phenomena  is  not  straightforward,  espe¬ 
cially  when  they  have  to  be  implemented  in  a  simulation 
tool.  The  time  constant  characterising  the  dynamic  behaviour 
of  the  battery  (e.g.  the  battery  voltage)  is  of  the  same  or¬ 
der  of  magnitude  as  the  mechanical  driving  time  constants. 
A  dynamic  battery  model  can  thus  lead  to  a  more  accurate 
simulation  result.  Different  models  have  already  been  pro¬ 
posed.  Fig.  6  illustrates  an  example  of  a  dynamic  battery 
model  for  a  lead  acid  battery  [28].  Fig.  7  shows  a  possible 
equivalent  circuit  for  a  NiCd  battery  [29].  Specific  battery 
management  and  test  infrastructure  has  recently  been  devel¬ 
oped  at  the  Vrije  Universiteit  Brussel  to  perform  dynamic 
battery  testing  in  a  well-controlled  environment.  With  the 
help  of  powerful  surface  fitting  algorithms,  dynamic  models 


Ubat 


Ri 


h 


U0 


as  well  as  a  new  type  of  SoC-indicator  are  under  develop¬ 
ment.  Up-to-now  a  maximum  open-loop  voltage  deviation 
in  time  is  implemented  in  the  software  model. 

With  some  types  of  batteries  the  electrolyte  can  start 
gassing  when  the  battery  cell  voltage  exceeds  a  certain  level, 
Ug,  and  the  charge  acceptance  goes  down.  For  lead  acid 
batteries  gassing  occurs  above  2.4  V  per  cell.  Whenever  the 
voltage  exceeds  this  value  a  constant  charge  acceptance  of 
80%  is  assumed  in  the  simulation. 

In  Eq.  (3.1 1)  a  square  root  is  taken.  If  the  argument  of  this 
square  root  is  negative,  it  means  that  the  battery  is  not  able  to 
deliver  the  required  power,  P^at-  The  maximum  power  that 
can  be  delivered  is  thus  given  by  Eq.  (3.12).  Remark  that 
according  to  Eq.  (3.11)  the  battery  is  able  to  accept  negative 
powers  of  any  amplitude  since  the  square  root  then  certainly 
is  a  real  value: 


Anax  — 


J] 2 

_o_ 

4Ri 


(3.12) 


Fig.  5.  Equivalent  scheme  of  a  simple  battery  model. 


When  the  battery  cell  voltage  exceeds  the  maximum  al¬ 
lowed  voltage,  (e.g.  during  regenerative  braking  on  a 
fully  charged  battery),  the  battery  power  is  limited  to  the 
corresponding  maximum  power,  Pmax(Umax)  (Eq.  (3.13)). 
And  vice  versa,  when  the  battery  cell  voltage  drops  under 
a  minimum  level,  (e.g.  during  strong  acceleration)  the 
corresponding  maximum  power,  Pmax(Umin),  is  calculated 
with  Eq.  (3.14).  Note  that  Pmax(Umax)  is  always  positive  and 
Emax(Umin)  always  negative.  The  possible  battery  power, 
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Fig.  6.  Dynamic  battery  model. 


Fig.  7.  Equivalent  circuit  for  the  NiCd  battery. 
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Table  2 


Speed  cycle:  simulation  vs.  measurements 


Measured 

Simulation 

ECE-15 

Ah  discharged 

1.52 

1.46 

1.53 

Ah  charged 

0.19 

0.18 

0.19 

Total  Ah 

1.33 

1.28 

1.34 

Recuperation  (%) 

12.35 

12.20 

12.71 

Ppos,  is  thus  limited  to  one  of  these  values  (Eq.  (3.15)). 
These  voltage  limits  Umax  and  Umm  are  required  in  real  life 
too.  Outside  these  limits  irreversible  electrochemical  pro¬ 
cess  will  damage  the  battery: 

tj  _  TSTrS  rrmax 

/W(£/max)  =  -  cell  (3.13) 


^max(f^min)  — 


U0  ~  Nrst/™f 


cell  TvT^Sj-^mm 


Ri 


Nr  £/; 


cell 


(3.14) 


E  pos  —  E 


max 


(f^max)  V  P max  (fAnin) 


(3.15) 


Furthermore  two  other  operating  limits  are  implemented  in 
the  simulation  programme:  the  maximum  battery  charge  cur¬ 
rent,  /^ax,  and  the  maximum  discharge  current,  /^x . 


3.3.  Calculation  results  of  the  battery  model 


To  illustrate  the  battery  model  some  simulations  as  well 
as  on-road  measurements  were  performed  for  a  passenger 
car  equipped  with  a  DC  separated  excited  motor  and  a  NiCd 
battery.  The  same  speed  cycle  as  the  one  driven  on  road  is 
simulated.  The  comparison  of  both  simulated  and  measured 
parameters  demonstrates  a  good  correlation  (see  Table  2). 
The  relative  error  is  less  than  5%. 

One  of  the  most  difficult  experiments  to  simulate,  while 
the  battery  model  plays  an  important  role,  is  an  acceleration 
test.  Contrary  to  a  comparison  based  on  a  pre-defined  speed 
cycle,  one  is  not  performing  a  straightforward  step-by-step 
calculation,  but  for  each  point  the  simulation  has  to  iterate 
towards  the  possible  working  point  [1].  A  little  error  in  the 
beginning  can,  due  to  integration,  result  in  a  large  deviation 
at  the  end  of  the  simulation. 

In  Fig.  8  one  can  find  the  measured  speed  compared 
against  the  simulation  results  when  the  acceleration  of  the 


Fig.  9.  Deceleration  simulation  vs.  measurement. 


vehicle  is  at  its  maximum.  The  boundaries  of  the  motor  are 
the  maximum  speed  and  torque.  This  motor  was  current  con¬ 
trolled.  A  current  limit  (as  a  function  of  the  revolutions  per 
minute)  is  introduced  too.  Simulated  values  are  marked  with 
‘s’  and  measured  values  with  Tn\  The  good  correlation  be¬ 
tween  the  measurement  and  the  simulation  demonstrates  the 
performance  of  the  iteration  algorithm.  An  average  devia¬ 
tion  of  2%  is  found. 

In  Fig.  9  one  can  see  the  deceleration  test.  During  this  test 
only  regenerative  braking  by  the  motor  was  performed,  with¬ 
out  using  any  mechanical  brakes.  The  graph  of  Fig.  10  com¬ 
pares  the  DC  motor  current  (Imot)  and  DC  voltage  (Umot) 
for  the  acceleration  and  deceleration  test.  The  little  devia¬ 
tion  in  current  can  be  explained  by  a  possible  minor  frontal 
or  back  wind  and  by  a  road  inclination  during  the  on-road 
measurement.  The  graph  of  Fig.  11  shows,  for  the  same  ac¬ 
celeration  test,  the  variations  in  current  (Ibat)  and  voltage 
(Ubat)  of  the  battery. 
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Fig.  10.  DC  motor  current  and  voltage  comparison. 
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Fig.  8.  Acceleration  simulation  vs.  measurement. 


Fig.  11.  Battery  voltage  and  current. 
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4.  Ultracapacitor 

4.1.  Ultracapacitor  description 

Super  or  ultracapacitors  (800-1500  F)  behave,  like  very 
high-power,  low-capacity  batteries  but  store  electric  en¬ 
ergy  by  accumulating  and  separating  opposite  charges 
physically,  as  opposed  to  batteries,  which  store  energy 
chemically  in  reversible  chemical  reactions.  One  key  as¬ 
pect  of  super-capacitors  is  that  they  demonstrate  excellent 
life  cycle.  When  fully  developed  for  vehicles,  they  could 
be  expected  to  last  as  long  as  the  car.  This  is  because  it 
is  possible  to  cycle  ultracapacitors  very  quickly  without 
having  the  large  decrease  in  life  cycle  that  most  chemical 
batteries  experience.  They  also  have  a  high  cycle  efficiency 
(up  to  90%  [30])  compared  with  chemical  batteries.  The 
primary  obstacles  with  ultracapacitors  at  this  time  are  their 
low  specific  energy,  which  are  in  the  range  5-10Wh/kg 
and  the  large  possible  voltage  dispersion  between  individ¬ 
ual  cells.  Ultracapacitors  also  have  the  unique  feature  that 
their  voltage  is  directly  proportional  to  their  state-of-charge. 
Therefore,  either  their  operating  range  must  be  limited  to 
high  state-of-charge  regions,  or  control  electronics  must 
be  used  to  compensate  for  the  widely  varying  voltage.  A 
voltage  drop  of  50%  can  be  allowed,  which  corresponds  to 
75%  of  the  energy  content.  The  fact  that  super-capacitors 
can  provide  high  power  for  accelerations  and  can  accept 
high  power  during  regenerative  braking  makes  them  ideally 
suited  for  the  load  levelling  required  in  electric  and  hybrid 
electric  vehicle.  Power  densities  of  2000-4000  W/kg  have 
been  demonstrated  by  ultracapacitors  in  the  laboratory  [31]. 

4.2.  Ultracapacitor  model 

The  model  of  an  ultracapacitor  can  be  developed  in  a 
similar  way  as  that  of  the  battery.  The  voltage,  Uc,  over 
the  capacity,  C,  has  the  same  roll  as  the  open-loop  voltage, 
t/0,  of  the  battery.  The  capacity  voltage,  Uc,  changes  with 
the  electrical  charge,  q.  The  electrical  charge  will  change  in 
function  of  the  cell  current,  /ceii ,  [32]: 

Uc  =  |  (4.1) 

4t+ 1  —  4t  fcell  T,  (4.2) 

The  equivalent  circuit  (Fig.  12)  is  also  similar  to  the  battery 
model. 

Several  capacitors  can  be  connected  in  series,  Nrs,  and 
in  parallel,  Nrp.  This  result  in  the  total  capacity  voltage, 
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Fig.  12.  Ultracapacitor  equivalent  circuit. 


f/cap,  and  current,  7cap  (Eqs.  (4.3)  and  (4.4)).  The  internal 
resistance,  R[,  is  dependent  on  the  current  and  temperature 
[33]: 


£/cap  =  Nr s(t/c  -  ft, /Ceii) 

(4.3) 

^cap  =  Nrp/ceii 

(4.4) 

Additionally,  the  self-discharge  of  the  capacitor  should  be 
taken  into  account  with  the  resistance,  Rv  [30] .  This  model 
does  not  take  into  account  the  possible  unbalance  between 
individual  cells.  This  model  can  be  used  as  a  black  box  model 
that  is  connected  to  the  propulsion  system  via  a  DC/DC 
converter. 


5.  Flywheel 

5.1.  Flywheel  description 

Flywheels  store  energy  mechanically  in  the  form  of  ki¬ 
netic  energy.  One  can  distinguish  mechanical  and  electrical 
flywheels.  The  mechanical  flywheel  is  connected  to  the  drive 
system  via  an  axle.  The  electrical  flywheel  takes  an  electri¬ 
cal  input  to  accelerate  its  rotor  by  using  a  built-in  motor,  and 
returns  the  energy  by  using  this  same  motor  as  a  generator. 

The  most  significant  factor  affecting  flywheel  design  is 
the  material  used  to  construct  the  flywheel  rim.  A  flywheel 
rim  needs  to  be  made  of  a  high  tensile  strength  to  density 
ratio  material  to  maximise  the  kinetic  energy  stored  (thanks 
to  a  high  rotational  speed)  while  minimising  the  chance 
of  failure.  Flywheels  show  technical  promise  as  a  potential 
load-levelling  device  for  hybrid  electric  vehicles;  however, 
they  are  still  under  development.  Flywheels  could  become  an 
excellent  high  power  density  storage  devices.  An  optimistic 
estimate  for  specific  power  is  ranging  between  2000  W/kg 
(short-term)  and  8000  W/kg  (long-term)  [34].  The  corre¬ 
sponding  specific  energy  is  ranging  from  4  to  50Wh/kg 
[35]. 

5.2.  Flywheel  model 

The  required  flywheel-convertor  power,  Preq,  will  be  di¬ 
vided  by  the  convertor  efficiency,  77,  including  convertor 
losses  and  all  flywheel  losses  except  the  friction  moment. 
This  required  power  is  compared  with  the  maximum  power, 
PY nax,  that  the  flywheel  can  deliver.  The  flywheel  current  ro¬ 
tational  velocity,  a>fw,  the  maximum  acceleration,  amSLX ,  the 
flywheel  inertia,  /,  and  the  maximum  torque,  Tmax,  define 
this  maximum  power: 

P in  ax  =  &>fw  max(rmax,  Jci  max  )  (5.1) 

The  flywheel  velocity  must  be  kept  above  a  certain  mini¬ 
mum  speed,  comin,  otherwise  the  flywheel  will  not  be  able 
to  deliver  power  and  below  a  maximum  velocity,  &>max,  to 
avoid  mechanical  damage. 
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The  maximum  velocity  defines  the  maximum  energy  con¬ 
tent  of  the  flywheel,  Em ax  (Eq.  (5.2)  which  corresponds  with 
a  SoE  equal  to  1).  To  avoid  confusion  with  the  battery  model, 
a  state  of  energy  is  used  in  the  flywheel  model  instead  of  a 
SoC.  Up  to  80%  of  this  energy  can  be  used  during  operation 
of  the  flywheel: 


E 


max 


,  )Z  T 

wma  xJ 

2.3600 


(5.2) 


Even  if  the  flywheel  is  not  delivering  power,  its  energy  con¬ 
tent,  £cons,  will  decrease  due  to  the  friction  moment,  7fr.  The 
energy  consumption,  £cons,  is  computed  as  function  of  the 
simulation  time  increment,  Ts ,  using  the  following  equation: 
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3600 
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(5.3) 


The  flywheel  energy  will  decrease  or  increase  as  function  of 
the  sign  of  the  flywheel  power,  Pfw  (Eq.  (5.4)): 

(5.4) 
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The  corresponding  flywheel  state  of  energy  and  angular  ve¬ 
locity  can  be  calculated  according  to  the  following  equa¬ 
tions: 

El 

SoE  =  — (5.5) 

Em  ax 


&>fw  — 


3600.2 

~J 


(5.6) 


where  is  the  flywheel  energy  at  instant  t  (Wh);  Eff~l  the 
previous  flywheel  energy  at  instant  t  (Wh);  &>fw  the  velocity 
(rad/s)  and  J  the  inertia  (kgm2). 


5.3.  Calculation  results  of  the  flywheel  model 


Fig.  13  shows  the  simulation  results  of  an  electric  drive 
with  a  flywheel  as  energy  source  delivering  the  traction 


power  for  driving  an  ECE-15  speed  cycle.  At  the  end  of  the 
cycle  the  flywheel  is  recharged  with  a  constant  power  of 
2500  W.  On  the  left  side  of  Fig.  13  one  can  find  settings  in¬ 
dicating  the  required  power  as  well  as  some  constant  model 
parameters  like  maximum  and  minimum  flywheel  speed  and 
its  maximum  deliverable  torque.  The  graph  in  the  middle  of 
Fig.  13  illustrates  the  flywheel  speed,  torque  and  power.  The 
flywheel  speed  changes  in  function  of  the  delivered  power. 

On  the  right  side  of  Fig.  13  some  other  constant  param¬ 
eters,  like  weight,  friction  moment,  converter  efficiency, 
flywheel  inertia  and  maximum  acceleration,  can  be  found. 
Finally,  on  the  left  side  the  output  parameters  can  be  found, 
namely  the  delivered  power  (inclusive  losses),  rotational 
speed,  efficiency  and  SoC  as  well  as  the  maximum  available 
power. 


6.  Engine-generator  group  or  APU 

6.1.  Engine 

The  engine  model  described  in  this  chapter  can  be  used 
in  a  classical  thermal  vehicle  as  main  motor  or  in  a  hybrid 
vehicle  as  part  of  the  APU.  It  is  mainly  based  on  look-up 
tables  of  engine  maps  and  calculates  the  emissions  and  fuel 
consumption. 

The  total  engine  torque,  7ice,  is  function  of  engine  in¬ 
ertia,  /,  vehicle  acceleration,  ra,  resistive  torque,  Tr,  and 
auxiliary  alternator  power,  Pa it,  as  shown  by  the  following 
equation: 

Tice  =  Tr  +  —  +  Ta  +  Ja  (6.1) 

Internal  combustion  engines  require  a  minimum  shaft  speed, 
below  which  they  will  stall.  This  is  implemented  in  the  clutch 
or  in  the  control  algorithm  of  the  hybrid  vehicle. 
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Fig.  13.  Simulation  result  of  a  fuel  cell  model. 
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If  the  engine  is  running,  with  the  engine  velocity  higher 
than  idle  speed  and  the  resistive  torque  different  from  zero, 
then  the  engine  is  in  normal  operating  mode  (see  further 
down).  If  not,  the  engine  can  be  in  off-mode  with  the  ‘key 
switch’  and  all  output  parameters  at  zero.  If  the  key  switch 
is  on,  with  the  engine  velocity  strictly  lower  than  idle  speed, 
then  the  engine  is  starting. 

The  user  can  select  a  start-up  penalty.  This  start-up  penalty 
has  to  be  measured  and  can  be  imported  as  a  cluster  of 
correction  factors  for  each  type  of  emission  and  for  the  fuel 
consumption.  This  can  be  done  for  different  engine  start-up 
temperatures.  Hence  the  additional  fuel  consumption  and 
emissions  can  be  taken  into  account  at  start-up.  If  the  user 
prefers  not  to  use  this  feature  or  the  data  is  not  available,  the 
idle  emissions  and  consumption  are  used  during  this  start-up 
phase. 

In  normal  operating  mode  the  maximum  torque  is  defined 
in  a  look-up  table  as  function  of  the  angular  velocity  and  the 
total  torque.  Once  the  possible  engine  torque  is  known  the 
emissions  (in  g/s)  and  fuel  consumption  (in  1/s)  can  be  cal¬ 
culated  on  the  basis  of  iso-emission  and  consumption  curves 
(in  g/kW  h).  Out  of  the  possible  speed  and  torque  values,  the 
emissions  are  calculated  using  bilinear  interpolation  of  the 
discretised  iso-emission  curves,  which  are  entered  in  several 
arrays.  These  curves  or  maps  describe  the  static  behaviour  of 
an  internal  combustion  engine.  This  works  when  transients 
are  negligible.  However,  if  transients  are  large  the  interpo¬ 
lated  results  may  be  wrong  with  several  orders  of  magnitude 
[36]. 

If  the  engine  torque  is  negative — the  engine  is  slightly 
braking  the  drive  train — a  small  (lower  than  idle  mode)  con¬ 
stant  value  or  zero  fuel  consumption  and  emission  is  as¬ 
sumed. 

The  engine  will  produce  additional  emissions  when  the 
engine  is  cold.  A  cold/hot  emission  ratio  can  be  introduced, 
which  is  a  function  of  ambient  temperature  and  trip  length 
[37-39].  The  user  can  select  if  the  motor  operates  in  a  warm 
mode  or  cold  mode.  In  the  case  of  a  cold  engine  an  additional 
penalty  is  implemented  based  on  an  empirical  model,  which 
describes  emissions  and  fuel  consumption  correction  factor 
as  function  of  the  ambient  temperature  and  fuel  type  [37]. 

The  mechanical  power,  Pice,  at  the  engine  shaft  (Eq.  (6.2)) 
as  well  as  the  power  of  the  consumed  fuel,  Pfuei  (Eq.  (6.5)) 
is  calculated.  The  last  one  is  a  function  of  the  fuel  density, 
D ,  and  fuel  specific  energy,  £fuei: 


Rice  —  ^ice^ice 

(6.2) 

fco  =  fco(tt>ice,  7ice,  t °) 

(6.3) 

fco 

Cons  = - 

(6.4) 

1000  D 

Rfuel  =  Cons  EfUQ\  x  10 

(6.5) 

where  7jce  is  the  total  engine  torque  (Nm);  Tr  the  resistive 
torque  (Nm);  Pa\t  the  alternator  power  (W);  ra  the  accel¬ 
eration  torque  (Nm);  <z>ice  the  angular  velocity  (rad/s);  Eice 


the  mechanical  shaft  power  (W);  t®  the  ambient  tempera¬ 
ture  (°C);  fco  the  fuel  consumption  (g/s);  D  the  fuel  density 
(kg/1);  Cons  the  fuel  consumption  (1/s);  Efuei  the  specific  en¬ 
ergy  (MJ/1)  and  Pfuei  the  fuel  power  (W). 

If  engine  data  corresponding  to  an  engine  with  a  certain 
required  size  is  not  available,  an  engine  can  be  selected  out 
of  the  database  and  its  characteristics  can  be  scaled  down  (or 
up)  by  the  help  of  a  scaling  factor  as  a  first  approximation. 

Notwithstanding  that  the  engine  model  considers  cold  op¬ 
eration  and  start-up  penalties,  it  is  still  a  rather  simple  model. 
Several  other  parameters,  not  taken  into  account  up-to-now, 
can  influence  the  fuel  consumption  and  emissions. 

The  model  does  not  consider  any  time  delay  associated 
with  fuel  transport  or  torque  production.  The  model  does 
not  simulate  cranking  or  associated  emissions.  No  transients 
of  any  kind  are  included.  Mean  value  torque  model  is  used, 
corresponding  to  the  average  torque  over  the  engine  cycle. 

Recent  investigations  carried  out  in  the  frame  of  the  Ger¬ 
man  and  Swiss  emission  factor  programme  have  shown  that 
it  is  necessary  to  take  into  account  the  influence  of  altitude 
when  determining  pollutant  emissions  for  the  vehicle  con¬ 
cepts  investigated  in  cases  where  a  major  proportion  of  the 
mileage  is  on  roads  at  a  high  elevation  [39]. 

In  the  case  of  conventional  spark-ignition  and  diesel  vehi¬ 
cles,  the  emission  behaviour  generally  deteriorates  within  a 
service  interval.  In  the  case  of  catalyst  vehicles  an  unavoid¬ 
able  deterioration  in  the  degree  of  conversion  by  the  catalyst 
(due  to  thermal  ageing  and  contamination)  leads  to  an  in¬ 
crease  in  emission  with  increasing  mileage.  Defects,  tuning 
errors  and  lack  of  maintenance  are  superimposed,  in  prac¬ 
tical  operation,  on  the  physically  determined  reduction  in 
the  degree  of  conversion  by  the  catalyst.  Consequently,  the 
implemented  look-up  tables  describe  the  engine  behaviour 
corresponding  to  the  age  of  the  characterised  engine. 

Hydrocarbon  emissions  from  motor  vehicles  are  caused 
by  two  major  sources:  exhaust  emissions  and  evaporative 
losses  through  the  vehicle’s  fuel  system  (storage  tank,  car¬ 
burettor  or  injection  system,  fuel  pipes).  Evaporative  emis¬ 
sions  occur  as  a  result  of  fuel  volatility  combined  with  the 
variation  of  the  ambient  temperature  during  a  24-h  period 
or  the  temperature  changes  of  the  vehicle’s  fuel  system  that 
occur  over  a  normal  driving  procedure. 

The  model  could  be  further  improved  by  taking  into 
account  certain  high-emission  events.  Such  high-emission 
events  occur  during  phases  of  extreme  high  acceleration 
and  during  gear  changes.  The  duration  of  such  an  event  is 
usually  only  a  few  seconds,  but  the  emissions  level  might 
reach  many  (even  up  to  20,000)  times  the  level  of  emissions 
during  normal  operation.  This  is  especially  true  for  modern 
gasoline  engines  with  closed-loop  catalytic  converters  [33]. 

6.2.  Generator 

The  purpose  of  the  generator  model  is  to  start  the  en¬ 
gine  and  to  control  the  velocity  of  the  engine  as  function 
of  a  certain  optimisation  criteria  in  a  series  hybrid  electric 
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propulsion  system.  The  input  parameter  of  the  generator  is 
the  required  DC  power  that  it  has  to  deliver  to  the  DC-bus.  If 
one  would  use  an  AC-generator  the  same  black  box  approach 
is  possible,  where  the  generator  and  the  required  rectifier  are 
gathered  in  one  model.  The  generator  model  differentiates 
three  operating  modes:  OFF,  starting  and  normal  operating. 

The  generator  model  will  switch  the  engine  OFF  when 
this  is  required  by  the  basic  power  control  algorithm.  This 
can  only  be  done  when  the  starting  process  is  completely 
finished  and  when  the  user  did  not  select  to  use  idle  mode 
(instead  of  shutting  down  the  engine). 

Once  the  engine  is  switched  off  and  needs  to  restart  again 
or  if  it  was  never  before  engaged,  it  will  first  pass  the  start-up 
phase  before  it  can  deliver  the  required  power.  The  corre¬ 
sponding  starting  will  take  a  certain  time,  selectable  by  the 
user.  This  starting  time,  Tstart,  will  define  a  required  acceler¬ 
ation,  restart?  as  function  of  the  generator  initial  speed,  coin{, 
and  the  required  start-up  velocity,  upstart  (Eq.  (6.6)): 


<2  start 


^  start  &4ni 


^start 


(6.6) 


The  generator  velocity  will  increase  proportional  to  this 
acceleration  as  function  of  the  simulation  time  increment 
(Eq.  (6.7)): 

(6.7) 


^gen  ^ start  -^s  ~~ I b  ^gen 


Hence  the  generator  velocity  is  defined  and  the  resistive  out¬ 
put  torque,  Tr,  can  be  calculated  as  function  of  the  required 
generator  DC  power,  Pdc  (Eq.  (6.8)): 

Pdc 


Tr  = 


1)00 


(6.8) 


gen 


where  oo^en  is  the  generator  velocity  (rad/s);  co^  the  pre¬ 
vious  generator  velocity  (rad/s);  r)  the  generator  converter 
efficiency  (%)  and  Ts  the  time  increment  (s). 

Due  to  its  inertia,  /,  accelerating  the  generator  will  intro¬ 
duce  an  additional  inertia  torque,  Td  (Eq.  (6.9)).  Probably 
the  engine  will  not  be  able  to  deliver  this  torque,  due  to  the 
fact  its  velocity  is,  e.g.  lower  than  idle  speed.  At  this  moment 
the  battery  will  deliver  the  inertia  power  during  start-up: 


Td  =  Ja 


(6.9) 


Once  the  generator  velocity  is  higher  than  the  start-up  ve¬ 
locity  and  the  generator  is  not  switched  off,  the  normal  op¬ 
erating  mode  will  be  active.  In  this  mode  the  DC  power 
will  be  limited  to  the  maximum  generator  power,  Pm ax.  The 
model  will  also  impose  a  maximum  power  variation,  dP/dt. 
In  this  way  it  is  possible  to  simulate  moderate  variations 
of  the  engine  required  power.  Consequently,  the  additional 
fuel  consumption  due  to  the  dynamic  operation  of  the  en¬ 
gine  is  prevented.  Similar  to  the  power  variation  limitation, 
the  model  will  impose  a  given  maximum  velocity  variation 
or  generator  acceleration,  amax.  Hence  the  additional  inertia 
torque,  ra,  due  to  too  fast  acceleration,  can  be  limited. 

In  the  series  hybrid  vehicle,  the  engine  velocity  is  inde¬ 
pendent  from  the  vehicle  velocity.  This  makes  it  possible 


to  control  the  generator  speed  in  such  a  way  to  operate  the 
engine  following  different  criteria  like: 

•  At  constant  speed  independent  from  the  required  power; 

•  At  a  velocity  corresponding  to  the  minimum  fuel  con¬ 
sumption  of  the  engine  as  a  function  of  the  required  gen¬ 
erator  power; 

•  At  a  velocity  corresponding  to  the  minimum  NO*  emis¬ 
sions  (or  other  pollutant)  of  the  engine  as  function  of  the 
required  generator  power; 

•  At  a  velocity  corresponding  to  a  rule  that  is  a  compromise 
between  fuel  consumption  and  emission  minimisation  of 
different  pollutants.  An  energy  optimisation  does  not  au¬ 
tomatically  imply  optimisation  of  the  emissions. 

These  criteria  can  be  implemented  with  the  help  of  a 
look-up  table  that  describes  the  relation  between  generator 
power  and  generator  velocity. 

The  efficiency  can  be  defined  with  the  help  of  statistical 
surface  fitting  models  based  on  measured  data  as  function 
of,  e.g.  the  generator  power  and  velocity.  Comparable  to 
the  start-up  mode  the  resistive  torque  Tr  is  calculated  using 
Eq.  (6.8). 

The  generator  mechanical  power  is  calculated  using  the 
following  equation: 

Egen  =  (Er  +  Ta)&>gen  (6.10) 

6.3.  Calculation  results  of  the  generator  engine  model 

For  the  illustration  of  the  generator-engine  model  the  fol¬ 
lowing  assumptions  are  made: 

•  The  vehicle  is  driving  the  ECE  cycle. 

•  The  APU  consists  of  a  generator  connected  to  a  diesel 
engine. 

•  The  maximum  acceleration  of  the  generator  is  set  to 
4  rad/s2. 

•  The  maximum  APU-power  deviation  is  set  to  lOkW/s. 

Fig.  14  represents  an  example  of  a  series  hybrid  electric 
vehicle  (SHEV)  in  which  the  APU  delivers  a  constant  power 
of  24  kW  to  the  DC-bus.  This  power  corresponds  to  the 
highest  fuel  efficiency. 

The  battery  is  recharged  with  a  power  corresponding  to 
the  APU  power  minus  the  required  driving  power.  During 
braking,  the  brake  power  will  charge  the  battery  too  (see 
green  line).  After  90  s  the  battery  power  reaches  its  maxi¬ 
mum  allowed  value  (see  Fig.  14).  The  iteration  algorithm  of 
the  software  tool  will  reduce  as  function  of  the  power  man¬ 
agement  strategy,  the  APU  power,  to  keep  the  battery  power 
within  limits. 

In  the  beginning  of  the  speed  cycle  (first  3  s)  the  engine 
start-up  process  can  be  noticed.  A  look  of  this  process  is 
illustrated  with  Fig.  15. 

Before  the  APU  can  deliver  power  the  engine  velocity  is 
brought  above  idle  speed  (first  4  s).  The  battery  delivers  the 
required  power  (battery  power  (green  line)  is  positive)  via 
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Fig.  14.  SHEV-ECE-maximum  regeneration  power. 
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the  charger  (this  is  the  converter  connected  to  the  genera¬ 
tor).  This  power  corresponds  to  the  model  with  the  inertia  of 
the  APU.  The  start-up  acceleration  is  defined  by  the  start-up 
time  and  the  required  idle  speed.  Once  the  engine  velocity 
exceeds  the  idle  speed,  the  engine  will  further  accelerate  to 
the  required  operating  setpoint.  The  generator  control  pro¬ 
cess  defines  this  acceleration.  After  17  s  the  APU  delivers  a 
constant  power. 


7.  Conclusions 

Since  years  ago  the  automotive  industry  and  several  re¬ 
search  institutes  have  developed  simulation  models  to  eval¬ 
uate  vehicle  performance,  fuel  consumption  and  emissions. 
Most  models  are  oriented  towards  conventional  vehicles  with 
internal  combustion  engines. 

The  VSP  developed  by  the  Vrije  Universiteit  Brussel  also 
allows  the  simulation  of  these  drive  trains,  however,  it  is 
especially  develop  to  simulate  battery,  hybrid  and  fuel  cell 
electric  vehicles.  A  great  attention  is  devoted  to  the  modular¬ 
ity  of  the  simulation  models,  especially  the  different  energy 
sources  used  in  these  electric  drive  trains.  The  software  tool 
has  a  flexible  database  structure,  integrated  in  the  compo¬ 
nent  models,  allowing  an  easy  implementation  of  different 
kind  of  component  description,  in  the  form  of  look-up  ta¬ 


ble,  maps,  theoretical  equations,  and  empirical  formula,  in 
function  of  the  available  data. 

This  paper  gives  an  overview  of  different  energy  sources 
applicable  in  electric  and  hybrid  vehicles.  Simulation  mod¬ 
els  for  fuel  cells,  batteries,  ultracapacitors,  flywheels  and 
engine-generator  units  were  described.  These  components 
were  described  in  three  stages:  first  their  functionality  and 
characteristics  are  described;  next  the  way  these  characteris¬ 
tics  can  be  implemented  in  a  simulation  model  is  discussed 
and  finally  these  principles  are  illustrated  with  the  help  of 
some  calculation  results. 

Because  of  nontrivial  difficulties  to  access  data,  the  mod¬ 
els  are  a  trade-off  between  complex  detailed  accurate  mod¬ 
els  and  models  based  on  easy-available  data.  Nevertheless, 
the  use  of  detailed  data  is  possible.  The  simulation  results 
illustrate  the  ability  of  the  VSP.  The  models  of  the  different 
energy  sources  are  developed  in  such  a  way  that  they  are 
easily  interchangeable  and  hence  allowing  a  flexible  eval¬ 
uation  of  different  drive  train  structures  and  their  internal 
power  flows. 

A  good  correlation  between  measurements  and  simulation 
results  is  found,  which  demonstrates  the  high  simulation 
accuracy. 

Up  to  now  the  software  was  especially  develop  for  the 
evaluation  of  electric  vehicles  as  well  as  hybrid  drive  trains 
and  their  complex  power  managements  strategies.  Further 
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work  is  being  done  to  extend  the  component  database  as 
well  as  to  further  validate  the  software  and  this  especially 
in  the  field  of  state  of  the  art  fuel  cell  systems. 
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